密度泛函理论研究MN@H_2O (M=Ga,Ge, In,Sn, Sb; N=M,Al)团簇的特性

采用密度泛函理论的B3LYP, B3P86, B1B95, P3PW91和PBE1PBE方法结合SDD, LANL2DZ和CEP-121G基组计算了d~(10)组态二聚物MN(M=Ga, Ge, In, Sn和Sb; N=M和Al)的几何结构.采用B3P86/SDD进一步研究了MN@H_2O团簇的几何结构及吸附能.结果表明,水分子结合在二聚物M_2上时,对二聚物影响较大,对水分子自身影响较小.将M_2中Ga, Ge, In, Sn或Sb替换一个原子为Al时,水分子在GeAl和SnAl上的吸附能变化较大,而在GaAl, InAl和SbAl上吸附能变化较小.另外, H_2O吸附在Ga, Ge, In, Sn和Sb上时,与吸附在Al上时,吸附能的变化不大.     

A DFT study on mechanisms of CO2 coupling with propargylic alcohols using alkali carbonates

The mechanisms of CO₂ coupling with the propargylic alcohol using alkali carbonates M₂CO₃ (M = Li, Na, K, Cs) have been investigated by means of density functional theory calculations. The calculations reveal that the target product tetronic acid (TA) is yielded through two stages: (a) the formation of the α-alkylidene cyclic carbonate (αACC) intermediate via Cs₂CO₃-mediated carboxylative cyclization of the propargylic alcohol with CO₂, and (b) the conversion of the αACC intermediate with Cs₂CO₃ to produce the cesium salt of the TA. Since the overall kinetic barriers for the two stages are comparable and affordable, the excellent chemoselectivity to the TA should be primarily originated from the high thermodynamic stability of the cesium salt of the TA. Moreover, relative to the TA, the possibility to yield the by-product acyclic carbonate can be excluded due to the both kinetics and thermodynamic inferiority. This result is different from the organic base-mediated reaction. Alternatively, our calculations predict that CsHCO₃ together generated with the cesium salt of the TA might also be an available mediating reagent for the incorporation of CO₂ with the propargylic alcohol. Compared to other alkali carbonates M₂CO₃ (M = Li, Na, K), the stronger basicity of Cs₂CO₃ and the lower ionic potential of cesium ion can raise the effective concentration of the αACC intermediate, and thus the conversion of the αACC intermediate into the cesium salt of the TA can be achieved with high yield.     

A new database and benchmark of the bond energies of noble-gas-containing molecules

We have developed a new database of structures and bond energies of 59 noble-gas-containing molecules. The structures were calculated by CCSD(T)/aug-cc-pVTZ methods and the bond energies were obtained using the CCSD(T)/complete basis set method. Many wavefunction-based and density functional theory methods have been benchmarked against the 59 accurate bond energies. Our results show that the MPW1B95, B2GP-PLYP, and DSD-BLYP functionals with the aug-cc-pVTZ basis set excel in predicting the bond energies of noble-gas molecules with mean unsigned errors (MUEs) of 2.0 to 2.1 kcal/mol. When combinations of Dunning's basis sets are used, the MPW1B95, B2GP-PLYP, DSD-BLYP, and BMK functionals give significantly lower MUEs of 1.6 to 1.9 kcal/mol. Doubly hybrid methods using B2GP-PLYP and DSD-BLYP functionals and MP2 calculation also provide satisfactory accuracy with MUEs of 1.4 to 1.5 kcal/mol. If the Ng bond energies and the total atomization energies of a group of 109 main-group molecules are considered at the same time, the MPW1B95/aug-cc-pVTZ single-level method (MUE = 2.7 kcal/mol) and the B2GP-PLYP and DSD-PLYP functionals with combinations of basis sets or using the doubly hybrid method (MUEs = 1.9-2.2 kcal/mol) give the overall best result.     

Designing spirobifullerene core based three-dimensional cross shape acceptor materials with promising photovoltaic properties for high-efficiency organic solar cells

The development of organic electron acceptor materials is one of the key factors for realizing high-performance organic solar cells (OSCs). Nonfullerene electron acceptors, compared to traditional fullerene acceptor materials, have gained much impetus owing to their better optoelectronic tunabilities and lower cost, as well as higher stability. Therefore, 5 three-dimensional (3D) cross-shaped acceptor materials having a spirobifullerene core flanked with 2,1,3-benzothiadiazole are designed from a recently synthesized highly efficient acceptor molecule SF(BR) ₄ and are investigated in detail with regard to their use as acceptor molecules in OSCs. The density functional theory (DFT) and time-dependent DFT (TDDFT) calculations have been performed for the estimation of frontier molecular orbital (FMO) analysis, density of states analysis, reorganization energies of electron and hole, dipole moment, open-circuit voltage, photo-physical characteristics, and transition density matrix analysis. In addition, the structure-property relationship is studied, and the influence of end-capped acceptor modifications on photovoltaic, photo-physical, and electronic properties of newly selected molecules ( H1-H5 ) is calculated and compared with reference ( R ) acceptor molecule SF(BR) ₄ . The structural tailoring at terminals was found to effectively tune the FMO band gap, energy levels, absorption spectra, open-circuit voltage, reorganization energy, and binding energy value in selected molecules H1 to H5 . The 3D cross-shaped molecules H1 to H5 suppress the intermolecular aggregation in PTB7-Th blend, which leads to high efficiency of acceptor material H1 to H5 in OSCs. Consequently, better optoelectronic properties are achieved from designed molecules H1 to H5 . It is proposed that the conceptualized molecules are superior than highly efficient spirobifullerene core-based SF(BR) ₄ acceptor molecules and, thus, are recommended to experiments for future developments of highly efficient solar cells.     

3d过渡金属掺杂M_2@Si_(20)(M=Sc-Zn)团簇的密度泛函理论研究

利用密度泛函理论(DFT)研究3d过渡金属掺杂硅团簇的几何结构和稳定性,计算了绝热电子亲和能和垂直电离能,内嵌双金属间距,自旋磁矩等.结果表明内嵌的Sc、Ti、V、Mn金属二聚体和十二面体硅笼构成了稳定的富勒烯结构,随着d电子数目的增加其内嵌的富勒烯构型有部分畸变,总体而言Si_(20)团簇掺杂双金属后稳定性得到了提高.     

Sliding mode control for a diffuse interface tumor growth model coupling a Cahn–Hilliard equation with a reaction–diffusion equation

We consider the sliding mode control (SMC) problem for a diffuse interface tumor growth model coupling a Cahn–Hilliard equation with a reaction–diffusion equation perturbed by a maximal monotone nonlinearity. We prove existence and regularity of strong solutions and, under further assumptions, a uniqueness result. Then, we show that the chosen SMC law forces the system to reach within finite time a sliding manifold that we chose in order that the tumor phase remains constant in time.     

基于不确定理论的铁路货运需求预测

合理预测铁路货运需求是铁路管理部门建设、运营等决策基础。为应对铁路货运需求的复杂变化,基于Pearson相关性分析方法筛选出铁路货运需求的七个具有关键影响的因素,并结合不确定理论建立不确定多元线性回归模型,相应的铁路货运预测结果由传统单一值变成可能的需求区间范围,更加符合处于不确定环境下的铁路货运需求实际情况。选取国家统计局2004~2016年相关数据进实证研究,并与回归模型以及BP模型的预测结果对比分析,实验表明不确定多元线性回归的预测结果更加精确。     

A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

A theoretical computational model of a plate in hypersonic flow

A theoretical model of an elastic panel in hypersonic flow is derived to be used for design and analysis. The nonlinear von Kármán plate equations are coupled with 1st order Piston Theory and linearized at the nonlinear steady-state deformation due to static pressure differential and thermal loads. Eigenvalue analysis is applied to determine the system’s stability, natural frequencies and mode shapes. Numerically time marching the equations provides transient response prediction which can be used to estimate limit cycle oscillation amplitude, frequency and time to onset. The model’s predictive capability is assessed by comparison to an experiment conducted at a free stream flow of Mach 6. Good agreement is shown between the theoretical and experimental natural frequencies and mode shapes of the fluid–structure system. Stability analysis is performed using linear and nonlinear methods to plot stability, flutter and buckling zones on a free stream static pressure vs temperature differential plane.